Scaling of solvation force in 2D Ising strip 
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The solvation force for the 2D Ising strip is calculated via exact diagonalization of the transfer 
matrix in two cases: the symmetric case corresponds to identical surface fields, and the antisym- 
metric case to exactly opposite surface fields. In the symmetric case the solvation force is always 
negative (attractive) while in the antisymmetric case the solvation force is positive (repulsive) at 
high temperatures and negative at low temperatures. It changes sign close to the critical wetting 
temperature characterizing the semi-infinite system. The properties of the solvation force are dis- 
cussed and the scaling function describing its dependence on temperature, surface field, and strip's 
width is proposed. 

PACS numbers: 05.50.+q, 68.08.Bc 



I. INTRODUCTION 

Properties of the solvation force in various condensed 
matter systems have been the subject of very intensive 
research during recent years 0, H S i, H, B 0, 0, S EH • 
Lattice models play special role among considered sys- 
tems. Although many important results were obtained in 
this field by different methods of taking into account fluc- 
tuations which determine the properties of the analyzed 
systems, the approach based on the exact evaluation of 
the partition function via the transfer matrix method still 
plays a distinguished role. Below we report our results 
on the properties of the solvation force for 2D Ising strip. 
They are obtained via exact diagonalization of the trans- 
fer matrix which is then followed by numerical solutions 
of equations for eigenvalues. 

We consider Ising model on a two-dimensional square 
lattice consisting of M rows and N columns. There is no 
bulk magnetic field acting on the system and there are 
two surface fields: hi acts on spins in the first row and 
hi acts on the M-th row. The Hamiltonian of this model 
takes the standard form 

N M-l 



7~L({ s i}) = X! s n,mS n ,m+l 

n—1 m—1 

N M 

n—1 m—1 

N N 



where s r , 



±1 with n 



1,. 



1. 



(1) 



,M 



denotes the spin located in the n-th column and m-th 
row, and J is the coupling constant. Periodic bound- 
ary conditions in the horizontal direction are imposed: 

s N+l,m = Sl,m- 



Our purpose is to determine the properties of the sol- 
vation force experienced by the system boundaries. In 
the following we consider two cases: the symmetric case 
(S) corresponds to hi = h 2 , and for the antisymmetric 
case (AS) one has hi = — h%. 

Although the strip of finite width experiences no 
phase transition, we shall often refer to two critical 
temperatures: the bulk critical temperature k B T c = 
2 J I In (l + \pl), and the wetting temperature T w which 
characterizes the critical wetting transition in a semi- 
infinite system with surface field h±. The wetting tem- 
perature fulfills the equation 



cosh 2^1 = cosh2if - smh2Ke~ 2K , 



(2) 



where K = J/ (/cb^w), Ki = hi/ (fegT w ), and k-g is the 
Boltzmann constant. For a small surface field hi this 
equation leads to 



j ' 1 



V2) In ( 1 



For hi > J one has T w 
transition. 



> 0. 

(3) 

and there is no wetting 



II. SOLVATION FORCE 

The dimensionless free energy per one column is de- 
fined as 



f{T,h u M) 



lim 



lim 



InZ 
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(4) 
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where Z is the canonical partition function for the Hamil- 
tonian {!]) which we evaluate using the exact transfer ma- 
trix method [HI, E2] • The free energy may be separated 
into three types of contributions 

/ (T, hi, M) = Mf h (T) + f s (T, hi) + f int (T, h u M), 

where /b is the bulk free energy per one spin [13[ and 
f s is the surface free energy; both /b and / s are M- 
independent. The remaining term /; nt describes the in- 
teraction between the system boundaries. It tends to 



2 




0.9 0.95 1 1.05 1.1 




-1 1 ' ' ' ' ' 1 ' 1 

0.7 0.75 0.8 0.85 0.9 0.95 1 1.05 1.1 
T/T c 



FIG. 1: Comparison of the solvation force / s ^ lv for T w = 0.8T C 
(hi 0.60J) (solid curve) and T w = (hi — J) (broken 
curve) for symmetric case hi = hi. Both curves correspond 
to the value M = 50. 

as M goes to infinity and from this term the solvation 
force originates. 

The solvation force is, in general, defined as the minus 
derivative of /i„t wrt M . In our case, because M takes 
only integer values, we use the definition 

/ solv (T, hi,M) = - [f int (T, hi, M + 1) - / int (T, hx,M)] , 

(6) 

which leads to the following expression 

/ so iv (T, hi, M) = f (T, hi,M) - f (T, hi, M + 1) - / b . 

(7) 

The solvation force for the 2D strip has been already 
analyzed in the T w = case [|[ corresponding to hi = J. 
Our analysis covers the whole spectrum T w > 0, i.e. 
hi < J. To calculate the solvation force we used the 
methods described in [ll], [l2|, [3 [H| • The complete anal- 
ysis (including also the inhomogeneous boundary fields) 
will be published elsewhere [l6| ; here we discuss only the 
main results. 

First we briefly discuss the symmetric case (S). The 
corresponding solvation force / s s olv is plotted in Fig|TJ for 
two cases: T w > (solid curve) and T w = (broken 
curve). The difference between these two functions is 
only of quantitative nature — decreasing the surface field 
hi, i.e. increasing the wetting temperature T w results 
in decreasing the absolute value of the solvation force. 
The minimum of the solvation force / s ^ lv is located at 

^min > ^c- 

For opposite surface fields (AS) the solvation force f^f v 
in the T w > case, i.e., hi < J, differs substantially from 
the T w = 0, i.e., hi — J case. FigJU presents a typical 
plot of /^f v as a function of T. For low temperatures 
this force is negative (attractive) and has a minimum at 
^rain < ^w- The solvation force is negative at wetting 
temperature and has a zero at T* > T w . Above T* the 
solvation force is positive (repulsive) and has a maximum 
at Tj^f x < T c (for M large enough). This remains in 
contrast with the T w — case in which the solvation 
force is positive for all temperatures. 



FIG. 2: The solvation force / ao f v as a function of temperature 
for T w = 0.8T C (i.e. hi « 0.60J) and M = 50 in the antisym- 
metric case hi — —h^. Note that the difference between the 
zero of the solvation force T* and the wetting temperature 
T w is not visible in this scale. 

Exactly at the bulk critical temperature T c the depen- 
dence of the solvation force on M has been found using 
conformal invariance (l7l . [l8| 

/ s s olv (T c , hi,M) = + O (1/Af 3 ) , (8) 

/sofv {T c , hi,M) — *^ + 0(l/M 3 ). (9) 

This result, in particular the universal values of the am- 
plitudes, has been reproduced by our analysis. 

In the rest of this paper we exclusively discuss the an- 
tisymmetric case (AS). 

First we concentrate on the temperature T* at which 
the solvation force becomes zero. We have found that 
for fixed hi > and for M — ► oo this temperature ap- 
proaches the wetting temperature T w exponentially 

Il-^L = A(h 1 )c- B ^ M , M^oo, (10) 

c 

where A (hi) and B (hi) are positive functions of the sur- 
face field hi We have found that for hi — > 0, i.e., 
T w — > T c one has 

lira A(hi) = 0. (11) 

hi— >0 

We note that this result is different from the corre- 
sponding result obtained within the mean field theory 
where T* is exponentially shifted below T w . It is also dif- 
ferent form the corresponding result obtained for the re- 
stricted solid-on-solid (RSOS) model, where T* is equal 
to T w [H. 

III. SCALING FUNCTION 

In this section we discuss the scaling function that de- 
scribes the behavior of the solvation force f^f v (T, hi , M) 
for large M and subcritical temperatures. The relevant 
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FIG. 3: Convergence of the scaling function. The function 
yM is evaluated for x = 1 and for several values of y and M 
presented on the plot. 



scaling function X(x) has already been proposed to de- 
scribe f^f v in the T w = case [f| 



f*f v (T, hx = J,M) 



M 2 



X 



M 



f" CO 



(12) 



The correlation length for 2D Ising model close to T c 



is r m = &t 



-+-i 



(T c - T) /T c > and 



where t 

= [41n(l + \/2)l • The scaling function Af(a;) can 
be obtained numerically from the transfer matrix spec- 
trum and some of its properties can be proved analyti- 
cally (see Eq.@ and [20j]), namely 



X(0) = 23tt/48, 
X(x) = 2tt 2 /x, for x 



oo. 



(13) 
(14) 



Here we would like to extend this result to hi < J, i.e. 
T w > case. 

For 2D Ising model the gap exponent Ai = 1/2 and 
for T < T c the following scaling behavior 



where 



AI 



y 



A hj 
k B T c t 1 / 2 



(15) 



(16) 



comes into play in the limit M — > oo with x and y fixed. 
This implies additionally the t — ► 0, hi — ► and T w — ► T c 

limits. The coefficient A = [(l + \/2) / In (l + \/2)] 1/2 
in Eq . (|16|) has been introduced such that the value y = 1 
corresponds to T = T w and then Eq.® is satisfied. For 
y < 1 Eq. (fl"5|) gives the solvation force f^f v below the 
wetting temperature and for y > 1 above T w . 

Eq. (|15|) may be rewritten in the form leading to the 
scaling function 



y(x,y) = lim y M (x,y), 



M 



(17) 



y(x,0.5) + 

y(x,i.5 x 

^(x,2.5) * 

y(x,3.5) 

X(X) . 



FIG. 4: Plots of the scaling function y. All plots were ob- 
tained from y M for M = 200. 
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FIG. 5: The scaling function y plotted for y = 0.5, 1.5, 2.5, 3.5 
and the scaling function X . All points were calculated for 
M = 200. 



where 



yM {x, y) 



(18) 



M 2 f t 



solv 



T c (l-^),fc B T cy (^) 1/2 ,M 



Fig[3] presents plots of y^ for fixed value of x = 1 and 
selected values of y € [0, 1.2] to exhibit the convergence 
of the series 3^m- Typically, for M large enough one has 



y M (x, y)=y (x, y) + + o (m- 2 ) 



(19) 



To estimate the values of the function y we evaluated J^v/ 
for M = 200. We note that a different way of obtaining 
the function y, based on the least squares method which 
allows to calculate the functions y and C in Eq. (fl"9]) . 
leads to similar results (the differences are not visible on 
the scale of our figures). The plots of function y are 
shown on Fig 01 
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Next we investigate the relation between the scaling 
functions y (x,y) and X (x). The function X (x) is cal- 
culated in the limit M — > oo, t — > with hi and Mt fixed. 
By applying this limit to Eqs (fT5|) and (fT2")l one gets 

«tf fx) = lim y (x,y) . (20) 

y— >oo 

The functions ^ (x) and 3^ (x, y) plotted for selected val- 
ues of y are presented on FigO Additionally we have 
found that X (x) — y (x, y) cx y~ 2 for large y. 

The scaling function y (x, y) changes its sign, see Figs 
IH [5] The zeros of the scaling function are denoted by 
y* (x), i.e., y(x,y*(x)) — 0. We have found that for 
large x the function y* (x) approaches 1 exponentially 
form above. This allows us to show that in the scaling 
limit M — > oo, hi — > 0, and M/if hxed one has 

where / (£) and g (£) are positive functions which can be 
determined via an implicit formula 

A C 1/2 = [/(C)] V V(/(C)/^)> (22) 
5(0 = Al-f(C). (23) 

We note that Eq. (f2~T]) is different form Eq. (fT0|) because 
different limiting procedures were applied in these two 
cases. 



IV. CONCLUSIONS 



In the symmetric case of identical surface fields the 
solvation force is negative (attractive) and has minimum 
at supercritical temperature. The solvation force calcu- 
lated for system parameters such that T w > differs only 
quantitatively from the one in case T w = 0, see FigJTJ 

In the antisymmetric case of opposite surface fields the 
solvation force is positive at high temperatures and neg- 
ative at low temperatures; it changes its sign at tempera- 
ture T* > T w , see FigJTJ In the case hi fixed and M — » oo 
the difference T* — T w approaches exponentially quickly 
in M. 

The scaling function y(x, y) was proposed to describe 
the behavior of the solvation force for T < T c in the limit 
hi — ► and M — * oo, see FigsSl [5] We checked that in 
the limit of high surface field (y — > oo in Eq. (fl~5]) ) this 
scaling function approached the scaling function describ- 
ing the scaling behavior of the for solvation force in the 
T w = case. In addition, the zeros of the scaling func- 
tion were investigated to find the formula for T* — T w in 
the hi — ► limit. 
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